Spin excitations in ~K 2 Fe i+x Se 5 : linear response approach 
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Abstract 

Using ab initio linear response techniques we calculate spin wave spectra in K2Fe4 +:E Se5, and find it to be in excellent 
agreement with a recent experiment. The spectrum can be alternatively described rather well by localized spin Hamiltonian 
restricted to first and second nearest neighbor couplings. We confirm that exchange coupling between nearest neighbor Fe 
magnetic moments is strongly anisotropic, and show directly that in the ideal system this anisotropy has itinerant nature which 
can be imitated by introducing higher order terms in effective localized spin Hamiltonian (biquadratic coupling). In the real 
system, structural relaxation provides an additional source of the exchange anisotropy of approximately the same magnitude. 
The dependence of spin wave spectra on filling of Fe vacancy sites is also discussed. 
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Because there seems to be a close, if poorly under- 
stood, connection between magnetic excitations and su- 
perconductivity in the recently discovered families of Fe- 
based superconductors, a great deal of attention has been 
paid to the elementary magnetic excitations in these 
systems p], @] . While the number of materials with super- 
conducting properties is rapidly growing, they all share in 
common a phase diagram with an antiferromagnetic re- 
gion immediately adjacent to the superconducting phase. 
AFM order in most parent compounds are characterized 
by rather steep spin waves (antiferromagnons) and rel- 
atively low Neel temperature (~140K) [J 0. On the 
other hand, iron selenides K2Fe4+ x Se5 discovered very 
recently, which we will call here a "245" system, is a 
particularly interesting case, because while its supercon- 
ducting state appears to be similar in many respects to 
some of the other families, has a much higher Neel tem- 
perature (Tjv~560-K") and very large magnetic moment 
(^3/is) [3j-|5[. The first neutron experiments, performed 
only recently @, show collective spin excitations some- 
what similar to excitations in the Fe-pnictides such as 
CaFe2As2- A highly debated topic of discussion is the 
anisotropy of the exchange coupling. In Ref. Q the au- 
thors fit a set of anisotropic Heisenberg model parame- 
ters to their data, and argued that at least three effective 
nearest neighbor (NN) exchange parameters are needed 
to fit observed spectra in a satisfactory manner. In an- 
other recent work [?} this exchange anisotropy has been 
fully attributed to biquadratic exchange, supporting the 
model introduced in Ref. [fj. However, neutron linear re- 
sponse experiments do not contain sufficient information 
to establish whether or not higher order terms of localized 
spin Hamiltonian or other interactions are responsible for 
this anisotropy. Information from the linear regime can- 
not unambiguously distinguish the dependence of linear 
response behavior on environment and terms originating 
from higher order. The unique answer can be obtained 
using non-collinear band structure calculations together 
with a linear response theory Q. 

Here we analyze spin excitations in the 245 system 
using non-collinear density functional theory in tan- 
dem with a linear response approach to determine the 



spin wave (SW) spectra of both the parent compound, 
K2Fe4Se5, and its modification by addition of Fe on the 
'vacancy' sites as described below. We accurately re- 
produce observed features in SW spectra of the parent 
compound from this first principles approach, without 
recourse to empirical or adjustable parameters. Thus we 
argue that the Heisenberg model extracted from this the- 
ory provides a good description of the static transverse 
susceptibility, and the exchange parameters derived 
ab initio from it, which have a physical meaning as the 
second variation in total energy with respect to spin ro- 
tations, are accurately given by the theory. We can also 
unambiguously address what happens when exchange pa- 
rameters, which can be computed to whatever range is 
desired, are truncated to a short range, e.g. just first and 
second neighbors (J1-J2 model). 

For computational convenience we adopt a multiple- 
scattering approach in the long wave approximation Q, 
which is reasonable provided the local Fe moment is suf- 
ficiently large. This is the case in K2Fe4+2;Se5, where 
the Fe local moment is both measured to be ^3/is, and 
predicted to be so in density-functional theory. 

The parent compound K 2 Fe 4 Se5 has a vacancy on the 
Fe sublattice for each formula unit, which if filled would 
have a composition K^FesSes and be structurally similar 
to the other Fe superconductors. The ordered magnetic 
phase of K 2 Fe4Se 5 is unusual: blocks of four Fe atoms are 
coupled ferromagnetically in a square; the squares are ar- 
ranged antiferromagnetically in a Neel like configuration, 
which we denote the "block Neel" magnetic structure. In 
what manner the vacancy Fe sites are filled in supercon- 
ducting material is a crucial issue, because the parent 
compound is predicted in density-functional theory to be 
an insulator with a bandgap of ~0.4 eV [l(J. 

We first consider magnetic excitations of the 
parent compound, with reciprocal lattice vectors 
Gi= (0.4 -0.2,0), G 2 =(0. 2,0.4,0), and G 3 =(0,0,o/c), in 
units of 2ir/a. This cell contains four formula units of 
K2Fe4Se5 with two units per plane. We used experi- 
mental lattice constants from Ref. pjj with a=3.914A, 
and c/ a— 3.587. Because the LDA is known to under- 
estimate Fe-chalcogen and Fe-pnictogen bond lengths in 
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FIG. 1: Spin-wave spectra in K2Fe4Ses calculated from the static susceptibility in the LDA. Panel (a) shows w(g), with the 
right frame for q along the F— Gi line. It corresponds to the experimental spectra shown in Ref. [(| (bottom panels of Fig. 3), 
whereas the left frame corresponds to the top panel of that Figure. The middle frame shows spectra along the F— (G1+G2) 
line. Panel (b) shows eigenvectors of the six optical modes at F, and their corresponding frequencies (meV). The perspective 
is from the 2 axis, so each point contains two Fe atoms. Arrows corresponding to spin rotations of each of the 16 Fe atoms. 
In some panels only 8 arrows are visible; this occurs when the rotations in the upper and lower planes are in phase. In other 
panels all 16 arrows are visible; for these modes the lower and upper planes are 180° out of phase. The right frame shows how 
the SW get modified when the exchange interactions are restricted to Ji and J2 only. 



TABLE I: Bond lengths (A) of K 2 Fe 4 Se 5 in the block 
Ne structures calculated from PBE and LDA functionals, 
bandgap (eV) and magnetic moment M (hb)- "LDA@PBE" 
refers to calculations using the Barth-Hedin (local) functional, 
but calculated at the PBE minimum-energy geometry. The 
shorter (longer) of the Fe-Fe bond length is the intrablock (in- 
terblock) distance. Mean PBE Fe-Fe and Fe-Se bond lengths 
are close to the experimental values. "25% Fe" refers to an 
LDAOPBE calculation with 25% doping of the Fe vacancy 
sites, as discussed in the text. 
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these families of superconductors, we use the PBE func- 
tional to relax the crystal structure in the block Neel 
magnetic state. The resulting relaxed structure, is seen 
in Table U to agree well with experimental data. On the 
other hand, the PBE functional tends to overestimate 
the magnetic moments and magnetic stabilization energy. 
Thus to calculate magnetic excitations, we use the LDA, 
but adopt the crystal structure predicted by the PBE 
functional. 

The full susceptibility x _1 (q, to) is prohibitively expen- 
sive to calculate for such a large system. Thus we adopt 
an approach which relies on the atomic sphere approxi- 
mation (ASA) following Ref. [L§] . To check ASA approx- 
imation to DFT, we carefully checked the energy bands, 
density of states, and magnetic moment in ASA against 
LDA results and found satisfactory agreement for all 
cases described here; for example the ASA bandgap and 
magnetic moment of the parent compound were found to 



be 0.40 eV and 2.85^ B , close to the LDA@PBE result of 
Table H 

Pair exchange parameters J,j were obtained from a 
Fourier transform of </(q); the latter was calculated on 
fine a g-mesh of 16x16x4 divisions in order to reliably 
extract Jjy to very distant neighbors. Resulting spin wave 
spectra for the parent compound K2Fe4Ses are shown in 
Fig. [1] Focusing first on the T point, six optical modes 
are seen: two degenerate pairs clustered around 100 meV, 
and two others near 230 meV. (We defer discussion of 
the mode near 40 meV for the moment.) The middle 
panel of Fig. Q] depicts their eigenvectors. The lowest 
pair and second pair are essentially the same in a given 
plane, but top and bottom planes rotate against each 
other for the low-energy pair, and in phase for the high 
energy pair. The first high-energy mode consists uniform 
rotations spins against parallel sheets on the xz axis; the 
second is similar but upper and lower planes are out of 
phase. 

Left and right frames can be compared directly with 
neutron data, top and bottom panels of Fig. 3 of Ref. [y] . 
All modes match experimental data very well, when 
broadening in the experimental spectrum is taken into 
account. The acoustic mode, with linear dispersion for 
small q of Fig.[T]has a maximum of 86 meV at g=(L/2)Gi, 
slightly larger than the measured value in Ref. [6j. The 
four low-energy optical modes are all "breathing" modes: 
spins rotate radially away from the center of their square. 
These four separate modes cannot be resolved in exper- 
iment; instead the measured data reflects a superposi- 
tion of these modes. It is nearly dispersionless along the 
T—Gi line, with wwllO meV. For comparison to exper- 
iment four optical modes were averaged and shown as a 
dashed line in Fig. [TJ and the result is seen to be the 
same as the neutron data, though at slightly higher en- 
ergy. Finally the two high-energy optical modes show a 
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TABLE II: Exchange energies Jk, in meV, for the Heisen- 
berg hamiltonian H=^2 Jij S» * Sj, of the parent compound 
K2Fe4Se5. Here k denotes the neighbor index. Jk and J' k re- 
fer to couplings between pairs of like spin and opposite spin, 
respectively. Also shown is the Neel temperature, estimated 
from the RPA [jj]. 
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weak q dependence, with w«240 meV. Essentially iden- 
tical behavior is observed in the neutron data along the 
(110) line, though u> is slightly smaller. The Neel tem- 
perature, Tjy, estimated within the RPA, is about 12% 
smaller than the observed value. This is in fact very good 
agreement with experiment, since the RPA is known to 
underestimate critical temperatures for a given Hamilto- 
nian, by roughly 10%. 

Finally, consider the mode appearing near 40 meV. If 
the two Fe planes were identical and uncoupled this mode 
would exactly coincide with the acoustic mode. Thus it 
is essentially the acoustic mode with q z — \/2-2ir/c, zone- 
folded to the T point, because the actual crystal contains 
two planes of Fe. The splitting of this and the true acous- 
tic mode are a result of coupling between the two planes 
through exchange parameters along J z coupling different 
planes. This shows that J z is small, but not so small that 
planes are decoupled. 

A great deal of speculation has arisen in the community 
about the size and range of parameters J. Much effort 
has been expended extracting effective exchange parame- 
ters Jij by fitting a Heisenberg model with a few (2 or 3) 
neighbors to the observed data, as was done in Refs. [6| 
and [7| . Since, as we have shown, SW can well described 
by the Heisenberg model derived ab initio, these results 
enable us to unambiguously address questions about the 
range and environment dependence of Jij. In particu- 
lar, it is debated whether a J\-Ji model is sufficient or 
more distant neighbors are required. Table |H] shows Jij 
for the first few neighbors, distinguishing intra-block and 
interblock couplings. As widely thought, the first two 
neighbors give by far the largest contribution to J. To 
quantify this effect, we recalculated the SW spectrum 
with J restricted to only nearest and second neighbor. 
This leaves four independent parameters since interblock 
J and intrablock J' are distinct (Table |TT| . The resulting 
SW spectra, shown in the right panel of Fig. [TJ look qual- 
itatively similar to the full calculation, except that some 
modes become degenerate, and low-energy modes soften 
by about 25%. The predicted Tjv drops by a comparable 



amount, from 494K to 399K. Perhaps most importantly, 
the "zone-folded" acoustic mode noted above becomes 
degenerate with the true acoustic mode, because all in- 
terplane interactions are now excluded by construction. 

Remarkably, corrections to the J\-Ji approximation 
receive almost no contribution from J3; rather, they orig- 
inate from J4-J7 (especially J|). The magnetic structure 
is stabilized by strong AFM first and second neighbor in- 
terblock coupling, and weakly so by nearest intrablock 
exchange. Frustration is present, since the AFM second 
neighbor intrablock coupling acts strongly to destabilize 
the magnetic order, but it is overcome because there are 
half as many such pairs as their interblock counterparts. 
The Table also shows what Jy would obtain if the struc- 
ture were constrained to an ideal geometry. Relaxation 
causes a small deformation of the Fe atoms in squares; 
there is also some dispersion in the Fe-Se bond length 
around the average value (Table |T| . This structural re- 
laxation is critically important for magnetic stability. As 
Table |TT] shows, the stability of the (intrablock) square 
is strongly enhanced by relaxation: the NN coupling is 
increased and frustration in the 2nd NN is strongly re- 
duced. Without relaxation some of SW frequencies be- 
come complex, indicating that the collinear state is not 
stable. 
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FIG. 2: Spin-wave spectra for K2Fe4.2sSe5 as described in the 
text. Actual cell is Ki6Fe34Se4o. 

As the filling of Fe vacancy sites very likely play an 
essential role in mediating superconductivity, how ele- 
mentary excitations are modified when these sites are 
partially occupied is a crucial question. We consider first 
25% occupation of Fe vacancies. To retain an overall 
zero-spin configuration it is necessary to double the size 
of the four formula unit cell: we construct a supercell 
with Gi = (0.3,0.1,0) and G 2 = (-0.1,0.3,0), and populate 
two of the eight vacancy sites with a pair of Fe^ and Fe^ 
atoms in such a way as to preserve a two-fold rotational 
symmetry. The lattice was relaxed using the PBE func- 
tional, with the resulting equilibrium bond lengths shown 
in Table Q] S Ws are shown in Fig. [2j The acoustic mode 
(now zone folded) is largely unchanged and uncoupled 
from the other modes: there is a linear dispersion at small 
q and a maximum near 80 meV for, e.g. q=(0.3,0.1,0) 
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and q=(0. 1,-0.2,0). The (zone-folded) low-energy opti- 
cal modes now show some new features. In addition to 
the modes in the 80-120 meV range, new and largely dis- 
persionless modes appear around 60 meV. The eigenfunc- 
tions of these modes are very complex, with strong ad- 
mixtures of both vacancy sites and host sites. The high- 
frequency modes are almost unchanged as "vacancy" Fe 
atoms participate very little in them. 

At 25% filling the magnetic structure is stable; indeed 
T N for K 2 Fe 4 Se 5 , predicted to be 494 K in the Tyab- 
likov approximation, is only slightly reduced, to 485K 
for K 2 Fe4.25Se5. We also studied a case with 50% filling. 
A collinear magnetic state was stabilized; however the 
linear response calculation of the SW spectra revealed 
imaginary frequencies, indicating that the collinear state 
is not stable. 

The exchange anisotropy, </(., is unusually large 
in this material. It is natural to associate Jk—J' k with 
biquadratic coupling, as was done in Refs.[14|. How- 
ever, lattice relaxations are very important in this case; 
compare the ideal to relaxed J in Table [TTJ Our calcula- 
tions predict that the intrablock NN Fe-Fe bond length 
is calculated to be smaller than the intrablock one by 
~8% (Table I); and perhaps more important, there is 
a significant dispersion in the Fe-Se bond lengths. The 
strictly electronic part of biquadratic contribution is still 
present and our calculations for the ideal 245 (no relax- 
ation) show that biquadratic coupling is about 30% of 
J\, This result was obtained by set of calculations of 
non-collinear configurations that enable us to go beyond 
linear response. We fit the total energy of the different 
non-collinear states to effective spin Hamiltonian with 
higher order terms and found a large cos 2 # contribution 
from first NN, while contribution from higher harmonics 
appeared to be smaller. Also biquadratic coupling from 
the second NN is smaller. Thus the environment depen- 



dence of J originates from two sources: purely electronic 
terms for a fixed lattice, and an "exchange-striction" con- 
tribution, originating from spin-configuration-dependent 
lattice relaxation. 

These non-collinear calculations point to the micro- 
scopic origin of the biquadratic interaction. It appears 
that major contribution to this 'biquadratic' energy is 
produced by the change of the amplitude of magnetic 
moments. M was found to vary by ^0.3/xb, even though 
its magnitude is large and thought to be well described 
by a rigid local-moment picture. The presence of unusu- 
ally large interaction between longitudinal and transverse 
degrees of freedom reflects the relatively large itinerant 
component of magnetic interactions in this system. It is 
notable that the interaction between these kinds of fluc- 
tuations lies outside a linear response description. The 
LDA is a mean-field theory and does not incorporate fluc- 
tuations directly. It can nevertheless capture some as- 
pect of fluctuations through the imposition of constraints 
(spin orientations). 

In conclusion, we demonstrated that ab-initio linear 
response method nearly perfectly describes the observed 
spin wave spectra in the complicated 245 system. The ex- 
perimental spectrum alternatively can be described rea- 
sonably well by localized spin Hamiltonian with two near- 
est neighbor couplings. We confirm the anisotropy of the 
exchange coupling between NN Fe magnetic moments 
and established that this anisotropy is associated with 
significant biquadratic coupling which appears in part be- 
cause of longitudinal fluctuations. Structural relaxation 
provides an additional source of the exchange anisotropy 
of approximately the same magnitude. 
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